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Decays of b hadrons into final states containing a D° meson and a muon are used to 
measure the bb production cross-section in proton-proton collisions at a centre-of-mass 
energy of 7 TeV at the LHC. In the pseudorapidity interval 2 < rj < 6 and integrated over 
all transverse momenta we find that the average cross-section to produce fe-flavoured or 
o3 ■ 6-flavoured hadrons is (75.3±5.4±13.0) /ib. 



Keywords: LHC, 6-hadron, cross-section, bottom production 
PACS: 14.65.Fy, 13.10.He, 13.75Cs, 13.85-t 



1 Authors are listed on the following pages. 



The LHCb Collaboration 



R. Aaij 23 , C. Abellan Beteta 35 '" 1 , B. Adeva 36 , M. Adinolfi 42 , C. Adrover 6 , A. Affolder 48 , 

M. Agari 10 , Z. Ajaltouni 5 , J. Albrecht 37 , F. Alessio 6 ' 37 , M. Alexander 47 , M. Alfonsi 18 , 

P. Alvarez Cartelle 36 , A. A. Alves Jr 22 , S. Amato 2 , Y. Amhis 38 , J. Amoraal 23 , J. Anderson 39 , 

R. Antunes Nobrega 22,fc , R. Appleby 50 , O. Aquines Gutierrez 10 , A. Arefyev 30 , L. Arrabito 53 , 

M. Artuso 52 , E. Aslanides 6 , G. Auriemma 22,/ , S. Bachmann 11 , Y. Bagaturia 11 , D.S. Bailey 50 , 

V. Balagura 30 ' 37 , W. Baldini 16 , G. Barber 49 , C. Barham 43 , R.J. Barlow 50 , S. Barsuk 7 , 

S. Basiladze 31 , A. Bates 47 , C. Bauer 10 , Th. Bauer 23 , A. Bay 38 , I. Bediaga 1 , T. Bellunato 20 '\ 

K. Belous 34 , I. Belyaev 23,30 , M. Benayoun 8 , G. Bencivenni 18 , R. Bernet 39 , R.P. Bernhard 39 , 

M.-O. Bettler 17 ' 37 , M. van Beuzekom 23 , J.H. Bibby 51 , S. Bifani 12 , A. Bizzeti 17 ' 3 , 

P.M. Bj0rnstad 50 , T. Blake 49 , F. Blanc 38 , C. Blanks 49 , J. Blouw 11 , S. Blusk 52 , A. Bobrov 33 , 

V. Bocci 22 , B. Bochin 29 , E. Bonaccorsi 37 , A. Bondar 33 , N. Bondar 29 ' 37 , W. Bonivento 15 , 

S. Borghi 47 , A. Borgia 52 , E. Bos 23 , T.J.V. Bowcock 48 , C. Bozzi 16 , T. Brambach 9 , 

J. van den Brand 24 , L. Brarda 37 , J. Bressieux 38 , S. Brisbane 51 , M. Britsch 10 , N.H. Brook 42 , 

H. Brown , S. Brusa 16 , A. Biichler-Germann 39 , A. Bursche 39 , J. Buytaert 37 , S. Cadeddu 15 , 
J.M. Caicedo Carvajal 37 , O. Callot 7 , M. Calvi 20 ^, M. Calvo Gomez 35 '" 1 , A. Camboni 35 , 

W. Cameron 49 , L. Camilleri 37 , P. Campana 18 , A. Carbone 14 , G. Carboni 21 '-? , R. Cardinale 19 '' 1 , 

A. Cardini 15 , J. Carroll 48 , L. Carson 36 , K. Carvalho Akiba 23 , G. Casse 48 , M. Cattaneo 37 , 

B. Chadaj 37 , M. Charles 51 , Ph. Charpentier 37 , J. Cheng 3 , N. Chiapolini 39 , A. Chlopik 27 , 
J. Christiansen 37 , P. Ciambrone 18 , X. Cid Vidal 36 , P.J. Clark 46 , P.E.L. Clarke 46 , 

M. Clemencic 37 , H.V. Cliff 43 , J. Closier 37 , C. Coca 28 , V. Coco 52 , J. Cogan 6 , P. Collins 37 , 

A. Comerma-Montells 35 , F. Constantin 28 , G. Conti 38 , A. Contu 51 , P. Cooke 48 , M. Coombes 42 , 

B. Corajod 37 , G. Corti 37 , G.A. Cowan 46 , R. Currie 46 , B. D'Almagne 7 , C. D'Ambrosio 37 , 

I. D'Antone 14 , W. Da Silva 8 , E. Dane' 18 , P. David 8 , I. De Bonis 4 , S. De Capua 21 J, 

M. De Cian 39 , F. De Lorenzi 12 , J.M. De Miranda 1 , L. De Paula 2 , P. De Simone 18 , D. Decamp 4 , 

G. Decreuse 37 , H. Degaudenzi 38 ' 37 , M. Deissenroth 11 , L. Del Buono 8 , C.J. Densham 45 , 

C. Deplano 15 , O. Deschamps 5 , F. Dettori 15 ' c , J. Dickens 43 , H. Dijkstra 37 , M. Dima 28 , 
S. Donleavy 48 , P. Dornan 49 , D. Dossett 44 , A. Dovbnya 40 , R. Dumps 37 , F. Dupertuis 38 , 
L. Dwyer 48 , R. Dzhelyadin 34 , C. Eames 49 , S. Easo 45 , U. Egede 49 , V. Egorychev 30 , 

S. Eidelman 33 , D. van Eijk 23 , F. Eisele 11 , S. Eisenhardt 46 , L. Eklund 47 , D.G. d'Enterria 35 '", 

D. Esperante Pereira 36 , L. Esteve 43 , E. Fanchini 20 ' 1 , C. Fiirber 11 , G. Fardell 46 , C. Farinelli 23 , 
S. Farry 12 , V. Fave 38 , G. Felici 18 , V. Fernandez Albor 36 , M. Ferro-Luzzi 37 , S. Filippov 32 , 

C. Fitzpatrick 46 , W. Flegel 37 , F. Fontanelli 19 '\ C. Forti 18 , R. Forty 37 , C. Fournier 37 , 
B. Franek 45 , M. Frank 37 , C. Frei 37 , M. Frosini 17 ' e , J.L. Fun gueirino Pazos 36 , S. Furcas 20 , 

A. Gallas Torreira 36 , D. Galli 14 - 6 , M. Gandelman 2 , P. Gandini 51 , Y. Gao 3 , J-C. Gamier 37 , 
L. Garrido 35 , D. Gascon 35 , C. Gaspar 37 , A. Gaspar De Valenzuela Cue 35 '" 1 , J. Gassner 39 , 
N. Gauvin 38 , P. Gavillet 37 , M. Gersabeck 37 , T. Gershon 44 , Ph. Ghez 4 , V. Gibson 43 , 
V.V. Gligorov 37 , C. Gobel 54 , D. Golubkov 30 , A. Golutvin 49 ' 30 ' 37 , A. Gomes 1 , G. Gong 3 , 

H. Gong 3 , H. Gordon 51 , M. Grabalosa Gandara 35 , V. Gracco 19,/l , R. Graciani Diaz 35 , 

L.A. Granado Cardoso 37 , E. Grauges 35 , G. Graziani 17 , A. Grecu 28 , S. Gregson 43 , G. Guerrer 1 , 

B. Gui 52 , E. Gushchin 32 , Yu. Guz 34 ' 37 , Z. Guzik 27 , T. Gys 37 , G. Haefeli 38 , S.C. Haines 43 , 

T. Hampson 42 , S. Hansmann-Menzemer 11 , R. Harji 49 , N. Harnew 51 , P.F. Harrison 44 , J. He 7 , 
K. Hennessy 48 , P. Henrard 5 , J. A. Hernando Morata 36 , E. van Herwijnen 37 , A. Hicheur 38 , 

E. Hicks 48 , H.J. Hilke 37 , W. Hofmann 10 , K. Holubyev 11 , P. Hopchev 4 , W. Hulsbergen 23 , 
P. Hunt 51 , T. Huse 48 , R.S. Huston 12 , D. Hutchcroft 48 , F. Iacoangeli 22 , V. Iakovenko 7 ' 41 , 

C. Iglesias Escudero 36 , C. Ilgner 9 , P. Ilten 12 , J. Imong 42 , R. Jacobsson 37 , M. Jahjah Hussein 5 , 
O. Jamet 37 , E. Jans 23 , F. Jansen 23 , P. Jaton 38 , B. Jean-Marie 7 , M. John 51 , D. Johnson 51 , 
C.R. Jones 43 , B. Jost 37 , F. Kapusta 8 , T.M. Karbach 9 , A. Kashchuk 29 , S. Katvars 43 , 

J. Keaveney 12 , U. Kerzel 43 , T. Ketel 24 , A. Keune 38 , S. Khalil 52 , B. Khanji 6 , Y.M. Kim 46 , 
M. Knecht 38 , S. Koblitz 37 , A. Konoplyannikov 30 , P. Koppenburg 23 , M. Korolev 31 , 
A. Kozlinskiy 23 , L. Kravchuk 32 , R. Kristic 37 , G. Krocker 11 , P. Krokovny 11 , F. Kruse 9 , 



ii 



K. Kruzelecki 37 , M. Kucharczyk 25 , I. Kudryashov 31 , S. Kukulak 25 , R. Kumar 14,37 , 

T. Kvaratskheliya 30 , V.N. La Thi 38 , D. Lacarrere 37 , G. Lafferty 50 , A. Lai 15 , R.W. Lambert 37 , 

G. Lanfranchi 18 , C. Langenbruch 11 , T. Latham 44 , R. Le Gac 6 , J. -P. Lees 4 , R. Lefevre 5 , 

A. Leflat 31 ' 37 , J. Lefrangois 7 , F. Lehner 39 , M. Lenzi 17 , 0. Leroy 6 , T. Lesiak 25 , L. Li 3 , Y.Y. Li 43 , 
L. Li Gioi 5 , J. Libby 51 , M. Lieng 9 , R. Lindner 37 , S. Lindsay 48 , C. Linn 11 , B. Liu 3 , G. Liu 37 , 

S. Lochner 10 , J.H. Lopes 2 , E. Lopez Asamar 35 , N. Lopez-March 38 , P. Loveridge 45 , J. Luisier 38 , 

B. M'charek 24 , F. Machefert 7 , I.V. Machikhiliyan 4 ' 30 , F. Maciuc 10 , 0. Maev 29 , J. Magnin 1 , 

A. Maier 37 , S. Malde 51 , R.M.D. Mamunur 37 , G. Manca 15 ' c ' 37 , G. Mancinelli 6 , N. Mangiafave 43 , 
U. Marconi 14 , R. Marki 38 , J. Marks 11 , G. Martellotti 22 , A. Martens 7 , L. Martin 51 , 

D. Martinez Santos 36 , A. Massafferri 1 , Z. Mathe 12 , C. Matteuzzi 20 , V. Matveev 34 , E. Maurice 6 , 

B. Maynard 52 , A. Mazurov 32 , G. McGregor 50 , R. McNulty 12 , C. Mclean 14 , M. Merk 23 , 
J. Merkel 9 , M. Merkin 31 , R. Messi 21j ', S. Miglioranzi 37 , M.-N. Minard 4 , G. Moine 37 , 

S. Monteil 5 , D. Moran 12 , J. Morant 37 , P. Morawski 25 , J.V. Morris 45 , J. Moscicki 37 , 

R. Mountain 52 , I. Mous 23 , F. Muheim 46 , K. Midler 39 , R. Muresan 38 , F. Murtas 18 , B. Muryn 26 , 

M. Musy 35 , J. Mylroie-Smith 48 , P. Naik 42 , T. Nakada 38 , R. Nandakumar 45 , J. Nardulli 45 , 

A. Nawrot 27 , M. Nedos 9 , M. Needham 38 , N. Neufeld 37 , P. Neustroev 29 , M. Nicol 7 , L. Nicolas 38 , 

S. Nies 9 , V. Niess 5 , N. Nikitin 31 , A. Noor 48 , A. Oblakowska-Mucha 26 , V. Obraztsov 34 , 

S. Oggero 23 , O. Okhrimenko 41 , R. 01deman 15 < c , M. Orlandea 28 , A. Ostankov 34 , B. Pal 52 , 

J. Palacios 39 , M. Palutan 18 , J. Panman 37 , A. Papadelis 23 , A. Papanestis 45 , M. Pappagallo 13 ' a , 

C. Parkes 47 ' 37 , C.J. Parkinson 49 , G. Passaleva 17 , G.D. Patel 48 , M. Patel 49 , S.K. Paterson 49 - 37 , 
G.N. Patrick 45 , C. Patrignani 19 '' 1 , E. Pauna 28 , C. Pauna (Chiojdeanu) 28 , 

C. Pavel (Nicorescu) 28 , A. Pazos Alvarez 36 , A. Pellegrino 23 , G. Penso 22 ' fc , M. Pepe Altarelli 37 , 
S. Perazzini 14 ''', D.L. Perego 20 '*, E. Perez Trigo 36 , A. Perez-Calero Yzquierdo 35 , P. Perret 5 , 
G. Pessina 20 , A. Petrella 16 ^ 37 , A. Petrolini 19 '\ E. Picatoste Olloqui 35 , B. Pie Vails 35 , 

D. Piedigrossi 37 , B. Pietrzyk 4 , D. Pinci 22 , S. Playfer 46 , M. Plo Casasus 36 , M. Poli-Lener 18 , 
G. Polok 25 , A. Poluektov 44 - 33 , E. Polycarpo 2 , D. Popov 10 , B. Popovici 28 , S. Poss 6 , 

C. Potterat 38 , A. Powell 51 , S. Pozzi 16 ' d , T. du Pree 23 , V. Pugatch 41 , A. Puig Navarro 35 , 
W. Qian 3,7 , J.H. Rademacker 42 , B. Rakotomiaramanana 38 , I. Raniuk 40 , G. Raven 24 , 

S. Redford 51 , W. Reece 49 , A.C. dos Reis 1 , S. Ricciardi 45 , J. Riera 35 '" 1 , K. Rinnert 48 , 

D. A. Roa Romero 5 , P. Robbe 7 ' 37 , E. Rodrigues 47 , F. Rodrigues 2 , C. Rodriguez Cobo 36 , 

P. Rodriguez Perez 36 , G.J. Rogers 43 , V. Romanovsky 34 , E. Rondan Sanabria , M. Rosello 35 '" 1 , 

J. Rouvinet 38 , L. Roy 37 , T. Ruf 37 , H. Ruiz 35 , C. Rummel 11 , V. Rusinov 30 , G. Sabatino 21 J, 

J.J. Saborido Silva 36 , N. Sagidova 29 , P. Sail 47 , B. Saitta 15 ' c , T. Sakhelashvili 39 , C. Salzmann 39 , 

A. Sambade Varela 37 , M. Sannino 19 ' ft , R. Santacesaria 22 , R. Santinelli 37 , E. Santovetti 21 '- 7 , 

M. Sapunov 6 , A. Sarti 18 , C. Satriano 22 '', A. Satta 21 , T. Savidge 49 , M. Savrie 1M , D. Savrina 30 , 

P. Schaack 49 , M. Schiller 11 , S. Schleich 9 , M. Schmelling 10 , B. Schmidt 37 , O. Schneider 38 , 

T. Schneider 37 , A. Schopper 37 , M.-H. Schune 7 , R. Schwemmer 37 , A. Sciubba 18,fc , M. Seco 36 , 

A. Semennikov 30 , K. Senderowska 26 , N. Serra 23 , J. Serrano 6 , B. Shao 3 , M. Shapkin 3 , 

I. Shapoval 40 ' 37 , P. Shatalov 30 , Y. Shcheglov 29 , T. Shears 48 , L. Shekhtman 33 , V. Shevchenko 30 , 

A. Shires 49 , S. Sigurdsson 43 , E. Simioni 24 , H.P. Skottowe 43 , T. Skwarnicki 52 , N. Smale 10,51 , 

A. Smith 37 , A.C. Smith 37 , N.A. Smith 48 , K. Sobczak 5 , F.J.P. Soler 47 , A. Solomin 42 , 

P. Somogy 37 , F. Soomro 49 , B. Souza De Paula 2 , B. Spaan 9 , A. Sparkes 46 , E. Spiridenkov 29 , 

P. Spradlin 51 , A. Srednicki 27 , F. Stagni 37 , S. Stahl 11 , S. Steiner 39 , O. Steinkamp 39 , 

O. Stenyakin 34 , S. Stoica 28 , S. Stone 52 , B. Storaci 23 , U. Straumann 39 , N. Styles 46 , 

M. Szczekowski 27 , P. Szczypka 38 , T. Szumlak 47,26 , S. T'Jampens , E. Tarkovskiy 30 , 

E. Teodorescu 28 , H. Terrier 23 , F. Teubert 37 , C. Thomas 51 ' 45 , E. Thomas 37 , J. van Tilburg 39 , 
V. Tisserand 4 , M. Tobin 39 , S. Topp-Joergensen 51 , M.T. Tran 38 , S. Traynor 12 , U. Trunk 10 , 
A. Tsaregorodtsev 6 , N. Tuning 23 , A. Ukleja 27 , O. Ullaland 37 , P. Urquijo 52 , U. Uwer 11 , 

V. Vagnoni 14 , G. Valenti 14 , R. Vazquez Gomez 35 , P. Vazquez Regueiro 36 , S. Vecchi 16 , 
J.J. Velthuis 42 , M. Veltri 17 '^, K. Vervink 37 , B. Viaud 7 , I. Videau 7 , X. Vilasis-Cardona 35 '" 1 , 
J. Visniakov 36 , A. Vollhardt 39 , D. Volyanskyy 39 , D. Voong 42 , A. Vorobyev 29 , An. Vorobyev 29 , 



iii 



H. Voss 10 , K. Wacker 9 , S. Wandernoth 11 , J. Wang 52 , D.R. Ward 43 , A.D. Webber 50 , 
D. Websdale 49 , M. Whitehead 44 , D. Wiedner 11 , L. Wiggers 23 , G. Wilkinson 51 , M.P. Williams 44 , 
M. Williams 49 , F.F. Wilson 45 , J. Wishahi 9 , M. Witek 25 , W. Witzeling 37 , M.L. Woodward 45 , 
S.A. Wotton 43 , K. Wyllie 37 , Y. Xie 46 , F. Xing 51 , Z. Yang 3 , G. Ybeles Smit 23 , R. Young 46 , 
O. Yushchenko 34 , M. Zeng 3 , L. Zhang 52 , Y. Zhang 3 , A. Zhelezov 11 , E. Zverev 31 . 

1 Centro Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 

2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 

3 Center for High Energy Physics, Tsinghua University, Beijing, China 
4 LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

5 Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
6 CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
7 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

8 LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 

9 Fakultdt Physik, Technische Universitdt Dortmund, Dortmund, Germany 

10 Max- Planck- Institut fur Kernphysik (MPIK), Heidelberg, Germany 

11 Physikalisches Institut, Ruprecht-Karls- Universitdt Heidelberg, Heidelberg, Germany 

12 School of Physics, University College Dublin, Dublin, Ireland 

13 Sezione INFN di Bari, Bari, Italy 

14 Sezione INFN di Bologna, Bologna, Italy 

15 Sezione INFN di Cagliari, Cagliari, Italy 

16 Sezione INFN di Ferrara, Ferrara, Italy 

17 Sezione INFN di Firenze, Firenze, Italy 

18 Laboratori Nazionali dell'INFN di Frascati, Frascati, Italy 

19 Sezione INFN di Genova, Genova, Italy 

20 Sezione INFN di Milano Bicocca, Milano, Italy 

21 Sezione INFN di Roma Tor Vergata, Roma, Italy 

22 Sezione INFN di Roma Sapienza, Roma, Italy 

23 Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands 

2A Nikhef National Institute for Subatomic Physics and Vrije Universiteit, Amsterdam, Netherlands 
25 Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland 
26 Faculty of Physics & Applied Computer Science, Cracow, Poland 
27 Soltan Institute for Nuclear Studies, Warsaw, Poland 

2S Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

29 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

30 'Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

31 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 

32 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 

33 Budker Institute of Nuclear Physics (BINP), Novosibirsk, Russia 

34 Institute for High Energy Physics (IHEP), Protvino, Russia 

35 Universitat de Barcelona, Barcelona, Spain 

36 Universidad de Santiago de Compostela, Santiago de Compostela, Spain 

37 European Organization for Nuclear Research ( CERN), Geneva, Switzerland 
3S Ecole Polytechnique Federate de Lausanne (EPFL), Lausanne, Switzerland 
39 Physik- Institut, Universitat Zurich, Zurich, Switzerland 

40 NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 

41 Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 

A2 H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 
^Department of Physics, University of Warwick, Coventry, United Kingdom 
i5 STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

46 School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 



iv 



47 School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 

48 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 

49 Imperial College London, London, United Kingdom 

50 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 

Department of Physics, University of Oxford, Oxford, United Kingdom 
52 Syracuse University, Syracuse, NY, United States 

53 CC-IN2P3, CNRS/IN2P3, Lyon- Villeurbanne, France, associated member 

5i Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 

a Universitd di Bari, Bari, Italy 

b Universitd di Bologna, Bologna, Italy 

c Universitd di Cagliari, Cagliari, Italy 

d Universitd di Ferrara, Ferrara, Italy 

e Universitd di Firenze, Firenze, Italy 

' Universitd di Urbino, Urbino, Italy 

9 Universitd di Modena e Reggio Emilia, Modena, Italy 

h Universitd di Genova, Genova, Italy 

% Universitd di Milano Bicocca, Milano, Italy 

3 Universitd di Roma Tor Vergata, Roma, Italy 

k Universitd di Roma La Sapienza, Roma, Italy 

1 Universitd della Basilicata, Potenza, Italy 

m LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain 

n Instituci6 Catalana de Recerca i Estudis Avangats (ICREA), Barcelona, Spain 



V 



1 Introduction 



Quantum Chromodynamics predicts the cross-section for the production of 6-flavoured 
hadrons in proton-proton collisions, for which higher order calculations are available [T]. 
The first data taken with the LHCb experiment at 7 TeV centre-of-mass energy allows 
this cross-section to be measured and compared to predictions. Knowledge of the b yield 
is also critical in ascertaining the sensitivity of experiments that aim to measure fun- 
damental parameters of interest involving, for example, CP violation. It is also useful 
for normalising backgrounds for measurements of higher mass objects that decay into bb, 
such as the Higgs boson. In this paper we present a measurement of the production cross- 
section for the average of 6-flavoured and 6-flavoured hadrons in proton-proton collisions 
at a centre-of-mass energy of 7 TeV in the pseudorapidity interval 2 < t] < 6, where 
rj = —In [tan(#/2)], and 9 is the angle of the weakly decaying b or b hadron with respect 
to the proton direction. We extrapolate this measurement to the entire rapidity interval. 
Our sensitivity extends over the entire range of transverse momentum of the fe-flavoured 
hadron. 

The LHCb detector [2] was constructed as a forward spectrometer primarily to measure 
CP violating and rare decays of hadrons containing b and c quarks. The detector elements 
are placed along the beam line of the LHC starting with the Vertex Locator (VELO), a 
silicon strip device that surrounds the proton-proton interaction region and is positioned 
8 mm from the beam during collisions. It provides precise locations for primary pp 
interaction vertices, the locations of decays of long lived hadrons, and contributes to 
the measurement of track momenta. Other detectors used to measure track momenta 
comprise a large area silicon strip detector located before a 3.7 Tm dipole magnet, and 
a combination of silicon strip detectors and straw drift chambers placed afterward. Two 
Ring Imaging Cherenkov (RICH) detectors are used to identify charged hadrons. Further 
downstream an Electromagnetic Calorimeter (ECAL) is used for photon detection and 
electron identification, followed by a Hadron Calorimeter (HCAL), and a system consisting 
of alternating layers of iron and chambers (MWPC and triple-GEM) that distinguishes 
muons from hadrons (MUON). The ECAL, MUON, and HCAL provide the capability of 
first-level hardware triggering. 

Two independent data samples, recorded at different times, are examined. For the 
earliest period of data taking the number of colliding bunches was sufficiently low that 
the high-level trigger could process all crossings and accept events when at least one track 
was reconstructed in either the VELO or the tracking stations. This data set, called 
"microbias", has an integrated luminosity, C, of 2.9 nb _1 . The second sample, referred to 
as "triggered", uses triggers designed to select a single muon. Here £ equals 12.2 nb -1 . 
These samples are analysed independently and the results subsequently combined. 

Most D° mesons are produced directly via pp — > ccX interactions, where X indicates 
any combination of final state particles. These particular D° mesons are denoted as 
"Prompt". D° mesons produced in pp — > bbX collisions where the 6-flavoured hadron 
decays into a final state containing a D° meson are called "Dfb". We use the decay 
channel b — >■ D°Xfi~V, as it has a large branching fraction of (6.84±0.35)% [3], and is 
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advantageous from the point of view of signal to background. Throughout this paper 
mention of a particular mode implies the inclusion of the charge conjugate mode as well. 

2 Analysis of D° K n + 

The Prompt and Dfb D° components can be separated statistically by examining the 
impact parameter (IP) with respect to the closest primary vertex, where IP is defined as 
the smallest distance between the D° reconstructed trajectory and the primary vertex H. 
We use the D° — > K~n + channel which has a branching fraction of (3.89±0.05)% [3]. 

The D° selection criteria are the same regardless of the trigger conditions. Both the 
kaon and pion candidates are associated with Cherenkov photons in the RICH system. 
The photon angles with respect to the track direction are examined and a likelihood 
formed for each particle hypothesis [2]. Candidates are identified as kaons or pions on 
the basis of this likelihood. We also require that the momentum transverse to the beam 
direction, p^, of both the kaon and pion be > 300 MeV, and that their scalar sum is >1400 
MeV. (We work in units with c=l.) Since real D° mesons travel before decaying, the kaon 
and pion tracks when followed backwards will most often not point to the closest primary 
vertex. We require that the \ 2 formed by using the hypothesis that the impact parameter 
is equal to zero, Xip> be > 9 for each track. They also must be consistent with coming 
from a common origin with vertex fit x 2 < 6. Finally, the D° candidate must be detached 
from the closest primary vertex. To implement this flight distance significance test we 
form a Xfs based on the hypothesis that the flight distance between the primary and D° 
vertices is zero, and require Xfs > 64. This set of requirements on the D° candidate is 
labeled "generic" . All of these requirements were selected by comparing sidebands of the 
invariant K~ir + mass distribution, representative of the background, with signal Monte 
Carlo simulation using PYTHIA 6.4 [5]. 

In order to ascertain the parameters characterizing the D° mass peak, a sample en- 
riched in Prompt D° mesons is selected. This is achieved by including two additional 
requirements: (1) the cosine of the angle between the D° candidate's momentum di- 
rection and the line from the K~ir + vertex to the primary vertex must be > 0.9999, 
and (2) the Xjp f° r the D° must be less than 25. The K~n + invariant mass distribu- 
tion after imposing all of these requirements is shown in Fig. [TJ The data are fit with a 
double- Gaussian signal function, with both Gaussians having the same mean, and a linear 
background. This signal shape is used in all subsequent fits. 

Selecting K~ir + candidates within ±20 MeV of the fitted D° mass peak and subtract- 
ing the background using invariant mass sidebands 35-75 MeV from the peak on both 
sides, we display the distribution of the natural logarithm of the D° candidate's IP in 
Fig. [2l Both Prompt and Dfb components are visible. The IP for the Prompt signal 
would be zero without the effects of resolution. The Prompt shape is described by a bi- 

Primary vertices are found using an iterative procedure based on the closest approach of tracks with 
each other. The resolutions of the resulting vertex positions depend on the number of tracks and are of 
the order of 70 /im along the beam direction and 10 ^m in each transverse coordinate, for 40 tracks. 
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Figure 1: K~tt + invariant mass for "Prompt" selection criteria in 2.9 nb" 1 . The curve shows 
a fit to a linear background (dashed) plus double-Gaussian signal function with parameters 
o"i=7.1±0.6 MeV, <72/<ti=1.7±0.1, and the fraction of the second Gaussian 0.40=L0.16. 
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Figure 2: Natural logarithm of the IP for D° mesons, with the IP in units of mm (points with 
error bars) for the 2.9 nb _1 microbias sample. Background has been subtracted using mass 
sidebands. The dashed curve shows the result of the fit to the Prompt component, the dotted 
line the Dfb component, and the histogram the sum of the two. 



furcated double- Gaussian function. The distribution for Dfb is widely spread as the finite 
b lifetime causes the D° meson not to point to the primary vertex; we use a Monte-Carlo 
simulated shape. The histogram in Fig. |2] shows the results of a fit to the two components, 
letting the parameters of the Prompt shape float; this shape is used in systematic studies. 
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3 Evaluation of the b 



D°Xfi-p yields 



3.1 Using microbias data 

To select the decay chain b — > D°Xfi~v, D° — > K~tt + and enrich our b sample, we match 
D° candidates with tracks identified as muons, by ensuring that they penetrate the iron 
of the MUON system and have minimum ionization in the calorimeters [2]. Right-sign 
(RS) combinations have the sign of the charge of the muon being the same as the charge 
of the kaon in the D° decay. Wrong-sign (WS) combinations have the signs of the charges 
of the kaon and the muon being opposite; they are highly suppressed in semileptonic b 
decay WS events are useful to estimate certain backgrounds. 

To find b candidates we select D° candidates using the generic criteria specified above, 
and add a track that is identified as a muon, has pt > 500 MeV, and has Xip > 4. 
The D° and muon candidates are required to form a common vertex with x 2 < 5, the 
D°n~ invariant mass must be between 3 and 5 GeV, and the cosine of the angle of the b 
pseudo-direction, formed from the D° and muon vector momentum sum with respect to 
the line between the D°fi~ vertex and the primary vertex, must be > 0.998. This angle 
cut is loose enough to have about 97% efficiency for b — > D°Xfi~V decays when taking 
into account the effect of the missing neutrino momentum. We measure t] using the line 
defined by connecting the primary event vertex and the vertex formed by the D° and the 
H~. Bins in rj are chosen to be larger than the resolution to obviate the need for any 
cross-feed corrections. Events are accepted in the interval 2 < rj < 6. 

The IP distributions of both RS and WS candidates, requiring that the K~n + invari- 
ant mass is within 20 MeV of the D° mass, are shown in Fig. [3j We perform an unbinned 
extended maximum likelihood fit to the two-dimensional distributions in K"ir + invariant 
mass over a region extending ±100 MeV from the D° mass peak, and ln(IP/mm). This 
fitting procedure allows us directly to determine the background shape from false combi- 
nations under the D° signal mass peak. The parameters of the Prompt IP distribution 
are found by applying the same criteria as for Fig. [3j but with the additional track failing 
the muon identification criteria. The Monte Carlo simulated shape is used for the Dfb 
component. 

The fit yields in the RS sample are 84.1±10.4 Dfb events, 16.3±5.4 Prompt events, 
and 14. Oil. 9 background. In the WS the corresponding numbers are 0.0±1.1 Dfb events, 
14.9±4.2 Prompt events, and 10.1±1.5 background. The Prompt yields are consistent 
between RS and WS as expected. 

The contribution of tracks misidentified as muons (fakes) in both the RS and WS sam- 
ples is evaluated by counting the number of tracks that satisfy all our criteria by forming 
a common vertex with a D° signal candidate, but do not satisfy our muon identification 
criteria. These tracks are categorized by their identity as electrons using ECAL, or pions, 
kaons or protons using the RICH. These samples are then multiplied by the relevant fake 
rates that were estimated from simulation and checked with data. The resulting ln(IP) 
distributions are examined, resulting in estimates of 2.2±0.4 RS Dfb fakes and 1.1±0.4 
WS Dfb fakes. The B(b D°Xt~V, r~ ji~uV) of (0.36±0.11)% is (5.3±1.6)% of the 
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Figure 3: Natural logarithm of the D° IP in the 2.9 nb 1 microbias sample for (a) right-sign 
and (b) wrong-sign -muon candidate combinations. The dotted curves show the D° sideband 
backgrounds, the thin solid curves the Prompt yields, the dashed curve the Dfb signal, and the 
thick solid curves the totals. 



semimuonic decay [3]. However, the relative efficiency to detect the resulting secondary 
muon is only 29% leading to a 1.5% subtraction. The lower efficiency is due to the lower 
secondary muon momentum from r decay and the finite r lifetime that causes some events 
to fail the vertex \ 2 requirement. Other sources of backgrounds from 6-hadron decays 
as evaluated by Monte Carlo simulation are small within our selection requirements, and 
predicted to be similar in size to the WS yields that are consistent with zero. 



3.2 Using muon triggered data 

The trigger imposes a cut of p^ > 1.3 GeV on muon candidates. The IP distributions 
for both RS and WS combinations are shown in Fig. HJ We find a total of 195.4±14.9 
RS Dfb, and 8.8±5.1 WS Dfb events. The Prompt contributions are determined to be 
9.3±4.8 RS with 5.3±3.0 WS. 

In order to extract the b cross-section from this data sample we have to make an 
additional correction for the overall //-dependent trigger efficiency. The Monte Carlo 
simulated efficiency is checked using data by studying J/ip — > decays in microbias 

events or those that triggered independently of the single muon trigger. The data show a 
somewhat larger relative efficiency than the simulation, from 2% at low 77 rising to 11% 
at high 77. We correct for this factor and use the 2% error determined on the correction, 
to account for its uncertainty, that we add to the statistical error of this sample. 

The IP distributions in each 77 bin in both trigger samples are fit independently to the 
same functions as described above to extract the 77-dependent event yields. The yields 
are listed in Table [TJ Muon fakes and the r~ contribution are subtracted in the same 
manner as in the microbias sample. In the triggered sample the hadron-to-muon fake rates 
are smaller as a result of the harder muon p^ cut imposed by the trigger of 1300 MeV 
rather than the 500 MeV used in analysing the microbias sample. The RS Dfb fakes total 
1.0±0.2 and the WS Dfb fakes total 0.6±0.2 events. A uniform 1.5% r~ subtraction is 
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Figure 4: Natural logarithm of the D° IP in the 12.2 nb 1 triggered sample for (a) right-sign 
and (b) wrong-sign -muon candidate combinations. The dotted curves show the D° sideband 
backgrounds, the thin solid curves the Prompt yields, the dashed curve the Dfb signal, and the 
thick solid curves the totals. 



done in each rj bin. 



Table 1: RS background subtracted event yields from data, and extracted cross-sections, com- 
pared with predictions from MCFM [7] and FONLL [8|. The systematic uncertainties in the 
normalisation of 17.3% are not included. The uncertainties on the FONLL prediction are ^33%, 
while those for MCFM are tf A %- 



Bin 


Event yields 




a(pp -> H b X) (rib) 








Microbias 


Trig. 


Microbias 


Trig. 


Average 


MCFM 


FONLL 


2 < 77 < 3 


16.7±4.5 


48.8±7.5 


27.2±7.3 


29.7±4.6 


29.0±3.9 


37.8 


28.9 


3 < rj < 4 


50.1±8.0 


111.4±11.0 


28.8±4.6 


28.8±2.8 


28.8±2.4 


27.1 


22.4 


4 < 7] < 5 


18.1±5.0 


30.2±6.0 


13.3±3.7 


11.7±2.3 


12.2±2.0 


16.7 


13.1 


5 < r] < 6 


4.7±2.8 


5.2±2.2 


6.5±3.6 


4.8±2.5 


5.3±2.0 


7.4 


5.9 


Sum 


89.6±10.8 


195.6±14.9 


75.9±10.0 


75.0±6.5 


75.3±5.4 


89.0 


70.2 



4 Luminosity determination and systematic uncer- 
tainties 

The luminosity was measured at specific periods during the data taking using both Van der 
Meer scans and the 'beam-profile' method [B]. Two Van der Meer scans were performed in 
a single fill. Analysis of these scans yielded consistent results for the absolute luminosity 
scale with a precision of around 10%, dominated by the uncertainty in the knowledge of 
the beam currents. In the second approach, six separate periods of stable running were 
chosen, and the beam-profiles measured using beam-gas and beam-beam interactions. 
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Using these results, correcting for crossing angle effects, and knowing the beam currents, 
we determine the luminosity in each period following the analysis procedure described in 
Ref. [6]. Consistent results were found for the absolute luminosity scale in each period, 
with a precision of 10%, again dominated by the beam current uncertainty. These results 
are in good agreement with those of the Van der Meer analysis. 

The knowledge of the absolute luminosity scale was used to calibrate the number 
of VELO tracks reconstructed using only the R sensors [2], which are found to have a 
stable response throughout the data-taking period. The integrated luminosities of the 
runs considered in this analysis were determined to be (2.85±0.29) and (12.2±1.2) nb _1 , 
respectively, for the microbias and triggered samples. 

The product of detector acceptance, tracking efficiencies and our analysis cuts, as 
estimated by Monte Carlo simulation, is about 8% for b hadrons produced in the region 2 < 
r] < 6. The systematic uncertainty on the tracking efficiency is evaluated by comparing 
the ratio of D° — > K~7t + 7t + tt~ to D° — > K~ir + events in data to the ratio in simulation. 
We find that the ratio of data to Monte Carlo efficiencies is 1.00±0.03 for tracks from D° 
decay, and use 3% as the uncertainty per track. For the higher momentum muon track 4% 
is used. The total tracking uncertainty then being fully correlated is taken as 10%, where 
this uncertainty is dominated by the size of the data sample. The kaon and pion RICH 
identification efficiencies are determined in each t] bin from a comparison of D° — > K~tt + 
yields evaluated both with and without kaon identification. An error of 1.5% is set on 
the particle identification efficiencies that is mostly due to the kaon, as the pion selection 
criteria are much looser. 

The efficiency of our muon selection criteria with respect to that obtained from the 
Monte Carlo simulation is evaluated as a function of momentum by detecting J /ip — >■ 
decays where one muon is identified by passing our muon identification criteria while the 
opposite-sign track must have been biased neither by the muon trigger, nor the muon 
identification criteria. Using the momentum weighted averages we find (data/MC) = 
(96.9^25)%>- We correct for the difference and assign a 2.5% error to our muon identifica- 
tion. 

Since the b — > D°Xfi~V detection efficiency changes linearly with by about a factor 
of four from zero to 12 GeV and then stays flat, the efficiency will not be estimated cor- 
rectly if the Monte Carlo generator does not accurately simulate the px distribution. We 
investigate this possible efficiency change by examining the difference between the mea- 
sured and simulated summed distribution of the D° plus muon. They are consistent, 
and an uncertainty of 3% is assigned as the systematic error from considerations of how 
large a difference the data allow. 

Because the detection efficiency is different for D° mesons that result from B~ — > 
D^°/j~V compared to those from other b decays (such as B — > D* + fi~V, B — > D**/_i~V, 
B s — > D**fi~V, or similarly from we include an uncertainty due to the precision of our 
knowledge of the branching fractions jl]. By varying these rates within their errors, we 
find an uncertainty of 4.4%. As discussed below, to translate our results on the yields into 
cross-section measurements we assume the fractions for fragmentation into the different 
6-hadron species as measured by LEP. Varying these values within their errors gives a 
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Table 2: Systematic uncertainties. 



Source 


Error (%) 


Source 


Error (%) 


Luminosity 


10.0 


Prompt & Dfb shapes 


1.4 


Tracking efficiency 


10.0 


B{D° ->K-ir + ) 


1.3 


B{b -> D°Xfj,-V) 


5.1 


D°fi~ vertex \ 2 cut 


1.2 


Assumed branching fractions 


4.4 


Kaon identification 


1.2 


LEP fragmentation fractions 


4.2 


Muon fakes 


1.0 


Generated b distribution 


3.0 


D° mass cut 


1.0 


Muon identification 


2.5 


D° vertex x 2 cut 


0.6 


Xip cut 


2.5 


D° flight distance cut 


0.4 


MC statistics 


1.5 


Pion identification 


0.3 


Total 




17.3% 





systematic uncertainty of 4.2%. 

The efficiency of the various selection criteria with respect to simulation has been 
evaluated by changing the cuts. The resulting changes of the yield are small. The D°p,~ 
vertex \ 2 cut efficiency was cross-checked comparing data and Monte Carlo using — > 
A/T - decays. All of the uncertainties considered are listed in Table [2J The total systematic 
uncertainty due to all sources added in quadrature is 17.3%. 



5 Cross-sections and comparison with theory 

The extracted cross-sections are listed in Table [TJ The //-dependent cross-section is shown 
in Fig. |5] for both data sets and the average. The agreement between the two data sets is 
excellent. 

We compare with two theories that predict b production cross-sections as a function 
of r\. MCFM [7] predicts the cross-section for bb quark production in next to leading 
order (NLO) using the MSTW8NL parton distribution function (PDF). The FONLL [S] 
prediction uses the CTEQ6.5 PDF, and improves the NLO result with the resummation 
of pt logarithms up to next-to-leading order. It also includes the b-quark fragmentation 
into hadrons. The measured yields are averaged over ^-flavoured and 5-flavoured hadrons, 
Hf,, in 7] intervals: 

# of detected _D°/i~ and D°fi + events 
<?{VV -> b ) = 2 £ x efficiency x B{b ->■ D°Xfj,-V)B (D° K~tv+) ' ( ' 

Averaging the cross-sections from both samples, and summing over t], we measure 

a(pp -> H b X) = (75.3 ± 5.4 ± 13.0) /xb (2) 

in the interval 2 < 77 < 6. The first error is statistical, the second systematic. The 
LEP fragmentation fractions are used for our central values [9]. Use of these fractions 
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Tl 

Figure 5: cr(pp — > H^X) as a function of r] for the microbias (x) and triggered (•) samples, 
shown displaced from the bin center and the average (+). The data are shown as points with 
error bars, the MCFM prediction as a dashed line, and the FONLL prediction as a thick solid line. 
The thin upper and lower lines indicate the theoretical uncertainties on the FONLL prediction. 
The systematic uncertainties in the data are not included. 

provides internal consistency to our results as B(b — > D°Xfi~V) was also measured at 
LEP. The measured value changes if the 6-hadron fractions differ. Fractions have also 
been measured at the Tevatron, albeit with large uncertainties [9]. The largest change 
with respect to LEP is that the 6-baryon percentage rises from (9.1±1.5)% to (21.4±6.8)%. 
If the Tevatron fractions are used, our result changes to (89.6±6.4±15.5) /A). 

6 Conclusions 

The cross-section to produce ^-flavoured hadrons is measured to be 

a(pp ->• H b X) = (75.3 ± 5.4 ± 13.0) /xb (3) 

in the pseudorapidity interval 2 < 7] < 6 over the entire range of pt assuming the LEP 
fractions for fragmentation into fe-flavoured hadrons. For extrapolation to the full r] region, 
theories predict factors of 3.73 (MCFM), and 3.61 (FONLL), while PYTHIA 6.4 gives 
3.77. Using a factor of 3.77 for our extrapolation, we find a total bb cross-section of 

a(pp -> bbX) = (284 ± 20 ± 49) /ib (4) 

based on the LEP fragmentation results; using the Tevatron fragmentation fractions the 
result increases by 19%. The quoted systematic uncertainty does not include any contri- 
bution relating to the extrapolation over the rj range where LHCb has no sensitivity. 

The production of 6-flavoured hadrons has been measured in pp collisions in 1.8 and 
1.96 TeV collisions at the Tevatron. The earlier measurements at 1.8 TeV appeared to be 
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higher than the NLO theoretical predictions [TO]. More recent measurements by the CDF 
collaboration at 1.96 TeV are consistent with the NLO theory [TT]. The history has been 
reviewed by Mangano [12]. Here, with a large energy increase to 7 TeV, we find that the 
measured cross-section is consistent with theoretical predictions, both in normalization 
and //-dependent shape. 
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